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Abstract 

The decay width of Z° to B c meson is evaluated at the next-to-leading order(NLO) accuracy 
in strong interaction. Numerical calculation shows that the NLO correction to this process is re- 
markable. The quantum chromodynamics(QCD) renormalization scale dependence of the results is 
obviously depressed, and hence the uncertainties lying in the leading order calculation are reduced. 
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I. INTRODUCTION 



B c meson has been attracting lots of attention in recent years. As the sole heavy quark 
meson that contains two different heavy flavors, its unique property attracts more and wide 



interests. Ever since its first discovery at the TEVATRON 



f ] , till now, various investigations 



on its production and decays have been carried out in aspects of theory 2|-|4| and experiment 



[i 



5, 6]. B c meson provides an excellent platform for testing the Standard Model(SM) and 
effective theories, e.g., to see whether non-relativistic QCD(NRQCD) [7j] is suitable for 
such system or not. In foreseeable near future, the B c physics study at the Large Hadron 
Collider (LHC) will tell us more about the nature of this special heavy bound system. 

Of the B_r meson production, apart from the direct ones, the indirect yields, like in top 
js] and Z° P -II] decays, are also important sources. The process of Z° decays to B c has 
an advantage of low background, but also with the disadvantage of low production rate, 
which prompts the LEP-I experiment unable to observe the B c signature [9|. In the future, 
if a high luminosity, I0 34 cm~ 2 s~ 1 or higher, electron-position collider, e.g. Internal Linear 
Collider (ILC) [12| . can set up, it would then bepossible to study the B c meson indirect 
production in Z° decays. As estimated by Ref. 13], there would be 10 9 ~ 10 10 Z° events 
produced each year at the ILC. Such kind of high luminosity collider, the so-called "Z 
factory" [14| , will provide new opportunities for both electroweak study and hadron physics. 

The indirect production of B c meson in Z° decays was evaluated by several groups at 
the leading order (LO) in strong interaction js-ll]. It is well-known that in charm- and 
bottom-quark energy regions, the higher order corrections of strong interaction are usually 
big, sometimes even huge. In order to make a more solid prediction on the B c production 
in Z° decays, and to depress the energy scale dependence lying at the LO calculation, an 
evaluation on the next-to-leading order(NLO) correction is necessary, which is the aim of 
this work. 

The paper is organized as follows: after the Introduction, in section II we repeat the 
leading order calculation on the Z° to B c decay width. In section III, the NLO virtual 
and real QCD corrections to Born level result are performed. In section IV, the numerical 
calculation for the process at NLO accuracy is done. The last section is remained for 
summary and conclusions. 



2 



B r 



(a) 



(b) 




(c) 




FIG. 1: The leading order Feynman diagrams for B c production in Z° decays. 



II. CALCULATION OF THE BORN LEVEL DECAY WIDTH 



At the leading order in a s , there are four Feynman Diagrams for B c production in Z° 
decays, the Z (j>i) — > B c (p Q ) + c(p 5 ) + b(p 6 ), as shown in FigJTJ The momentum of each 



particle is assigned as: p x = p z , Pz = pi, Pi = Pc, P5 = Pc, Pe = Pb, Po = P3 + Pi, P3 



Pi- 



For b and c quark hadronization to B c meson, we employ the following commonly used 
projection operator 



v(p 3 )u(p 4 



1 



--^siPo + m b + m c ) x 



(1) 



2x/2 \^/(rn b + rn c )/2 —'J \^N ( 

Here, l c stands for the unit color matrix, and N c = 3 for QCD. The nonperturbative 
parameter ipB c {fy is the Schrodinger wave function at the origin of the bc bound states. In 
our calculation, the non-relativistic relation rriB c = m b + m c is also adopted. 

The LO amplitudes for B c production can then be readily obtained with above prepara- 
tions. They are: 



M a 



7ia s gij Bc (0)C F 5 jk 



^(pi)[(4sin 2 g w -3) + 3 75 
(^i- ^ 3 - m b )(p A + p 5 ) 2 



ilM + m Bc )^v(p 5 ) , (2) 



7ia s gip Bc {0)C F 6 jk _ . /(pi) (8 sin 9 W - 3) + 375 
= fi /q «(Pe)7 M 175 (A + m B c J— ( — - — wrz -, ^7^(Ps) , (3) 



7^oi s gip Bc (0)C F 8 jk \ /(Pi) (8sm W - 3) + 375 

6v3mB c cos0 w p (^1- p 5 - m c )(p 3 + p G ) 2 



and 



= « /o a u (Ps)—( , ylTu u \~ 7 *75 W) + m Bc 7 M w (p 5 ) . 5 

6y/dm Bc cos 9 W (p 4 +p 5 ) 2 {pi- p 6 - m b ) 

Here, j, k are color indices, Cp = 4/3 belongs to the SU(3) color structure. 9w is the 
Weinberg angle with the numerical value sin 2 6w — 0.23. 

The Born amplitude of the processes shown in FigJT] is then M. Born = M. a + -Mb + -M c + 
M-d-, and subsequently, the decay width at LO reads: 

dr Born = -Li V |AW„| 2 dPS 3 . (6) 

IVTLz 3 Z — ' 

Here, symbolizes the sum over polarizations and colors of the initial and final particles, 
I comes from spin average of initial Z° meson, dPS3 stands for the integrals of three-body 
phase space, whose concrete form can be expressed as: 

dPS '=3^4^ dS " < 7) 

where Si = (po+Ps) 2 = {.Pi—Ps) 2 and s 2 = (ps+Pe) 2 = (pi—Po) 2 are Mandelstam variables. 
The upper and lower bounds of the above integration are 



mm 
in 



yj /VI, s 2, m|J • f[s 2 , m 2 c , m 2 } + \m\ -s 2 - {m b + m c ) 2 }(s 2 + m 2 c - m 2 b ) 
+m 2 Bc + ml , (8) 



\j f[m 2 z , s 2} mfj • f[s 2 , m 2 , m 2 } - [m| -s 2 - (m b + m c ) 2 ](s 2 + m 2 - m 2 b ) 

2i 2 

+m 2 + m 2 c (9) 



and 



[m z ~ (m b + m c )} 2 , s™ n = (m c + m b ) 2 (10) 
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with 



f[x,y,z] = [x -y - zf - Ayz . 



(11) 



III. THE NEXT-TO-LEADING ORDER CORRECTIONS 



At the next-to-leading order, the Z° boson decay to B c includes the virtual and real 
QCD corrections to the leading order process. For the two kinds of vertices, Z — >■ 66 and 
Z — > cc, we need only to consider one of them, e.g. as shown in FigsEHHJ since they are 
similar. For the virtual corrections, the decay width at the NLO can be formulated as 



dT Virtua i = —— - y^ 2Re(M* Born M V irtuai)dPS 3 . 
2mz 3 



(12) 



In virtual corrections, the ultraviolet (UV) and infrared(IR) divergences exist universally. 
We use the dimensional regularization scheme to regularize both UV and IR divergences, 
similar as performed in Ref. 15j, and use the relative velocity v to regularize the Coulomb 



divergence [16|. According to the power counting rule, the UV divergences exist merely 
in self-energy and triangle diagrams, which can be canceled by counter terms. The renor- 
malization constants include Z2, Z%, Z m , and Z g , corresponding to quark field, gluon field, 
quark mass, and strong coupling constant a s , respectively. Here, in our calculation the Z g 
is defined in the modified-minimal-subtraction (MS) scheme, while for the other three the 
on-shell (OS) scheme is adopted, which tells 
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(13) 



Here, the mass m in 5Z° S and bZ° s represents mj or m c ; (3 = (11/3)Ca — (4/3)T/n/ is 
the one-loop coefficient of the QCD beta function; rif = 5 is the number of active quarks in 
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FIG. 2: The self-energy diagrams in virtual corrections. 

our calculation; and (3' = (11/3)Ca — (4/3) T/ni/ with nif = 3 being the number of light- 
quark flavors; Ca = 3 and Tp = 1/2 attribute to the SU(3) group; /i is the renormalization 
scale. Note, since the terms related to <5Zg cancel with each other, the full NLO result is 
independent of the renormalization scheme of the gluon field. 

In dimensional regularization scheme, 75 is an object hard to handle; especially in the 
process that contains the vector- axial current, things become more complicated. In our 
work, we adopt the scheme provided in Ref.[r7|, where the following rules must be obeyed: 

I. the cyclicity is forbidden in traces involving odd number of 75. 

II. For the certain diagrams that contribute to a process, we must write down the am- 
plitudes starting at the same vertex, named the reading point. 

III. As a special case of rule II, if the anomalous axial current exists, the reading point 
of the anomalous diagrams must be the axial vector vertex, in order to guarantee the con- 
servation of the vector current. 

By utilizing this rule in our process, the two anomalous diagrams denoted as TriangleNll 
and TriangleN12 in Figj3]are calculated, and the UV divergences in these two diagrams are 
canceled by each other. To deal with the 75S except for what in anomalous diagrams, the 
cyclicity is employed to move the 75s together and then are contracted by 7! = 1. Hence, if 
a trace contains even number of 75, there will be no 75 left. Otherwise, after the contraction 
of odd number of 75, one remains. 
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FIG. 3: The triangle diagrams in virtual corrections. 



In the virtual correction, IR divergences remain in the triangle and box diagrams. Of 
all the triangle diagrams, only two have IR divergences, which are denoted by TriangleN3 
and TriangleN8 in Figj3) Of the diagrams in FigJU BoxN3 has no IR divergence, BoxN6 has 
merely the Coulomb singularity, PentagonN8 has both a Coulomb singularity and ordinary 
IR divergence, and the remaining other diagrams have only the ordinary IR divergences. 
We find that the combinations of BoxNl + BoxN4, BoxN5 + PentagonN9 + TriangleN8, and 
BoxN7 + TriangleN3 are IR finite, while the remaining IR singularities in BoxN2 and BoxN8 
are canceled by the corresponding parts in real corrections. The Coulomb singularities 
existing in BoxN6 and PentagonN8 can be regularized by the relative velocity v. The ^ 
terms are renormalized by the counter terms of external quarks which form the B c , while 
the - term will be mapped onto the wave function of the concerned heavy meson. In the 
end, the IR and Coulomb divergences in virtual corrections can be expressed as 



dT; 



\I R,Coulomb 
virtual 



dr 



8a s 7r 2 1 (m| + 2mt,m c — 2pi ■ po)x s lnx s 1 



(14) 
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FIG. 4: The box and pentagon diagrams in virtual corrections. 

. , , 1- a/ 'l-4m b m c / (4m b m c +m 2 z -2pi-po) TT , . i 1 • 

with pi = pz, Po = Pb c and x s = v , = Here, m this work - m 

l+-y / l-4m i ,m c /(4m i ,m c +m|-2pi-po) € 

fact represents - — Je + ln(47r/^ 2 ). 

Of the concerned process, there are 12 different Feynman diagrams in real correction, as 
shown in FigJSl Among them, RealN2, RealN3, RealN8, and RealN9 are IR-finite, meanwhile 
the combinations of RealNl + RealN5 and RealNIO + RealN12 exhibit no IR singularities as 
well, due to the reason of gluon connecting to either b or c quark of the final B c meson. The 
remaining diagrams, RealN4, RealN6, RealN7, and RealNll are not IR singularity free. To 
regularize the IR divergence, we enforce a cut on the gluon momentum, the p^. The gluon 
with energy p® < 5 is considered to be soft, while pj > 5 is thought to be hard. The 5 is a 
small quantity with energy-momentum unit. In this way, the IR term of the decay width 
can then be written as: 
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FIG. 5: The real correction Feynman diagrams that contribute to the production of B c . 



where dPS4 is the four-body phase space integrants for real correction. Under the condition 
of p° < S, in the Eikonal approximation we obtain 

d 3 p 7 



dPS 



4 \soft~- 



dPS, 



o \P° 7 <5 



(16) 



'(27r) 3 2p° 

In the small 8 limit, the IR divergent terms in real correction can therefore be expressed as 



dr 



IR 

Real 



dr 



Born ' 



3tt 



- - ln(<5 2 
e 



(m| + 2m b m c - 2p x ■ p )x s lnx s 



+ finite terms > . (17) 



mfem c (l — x 2 s ) 

Here, the ln(<5 2 ) involved terms will be canceled by the 5-dependent terms in the hard sector 
of real corrections. Comparing ( TT7|) with (fl4|) . it is obvious that the IR divergent terms 
in real and virtual corrections cancel with each other. In the case of hard gluons in real 
correction, the decay width reads 



j-phard 
Ui Real 



1 1 



Real\ 2 CIPS4 \hard 



2mz 3 

In this case, the phase space dPS 4 \hard can be expressed as 

yftsy + 



(18) 



dPS 



4 \hard~ 



Asyml 



(47T) f 



dp ° 



PQ- 



1 p2ir 

d cos 0c / d(f) c 
i Jo 



fP7- rv+ fP7+ ry+ 1 

/ dp 7 ° / dy + / dp 7 ° dy\ (19) 

■/<5 J Pr ° i (m " +mc) I 



with 



p = m b + m c , (20) 

P0° + = ^ , (21) 

Pr° = S ' 2V ~ SP °° > ( 22 ) 

o _ g ~ 2yW 

~ 77^= — " > ( 2 d) 



2v ^ _ 2 p o - 2^ 

V- = \\W*- Po° - Pr ) 2 - l^l 2 - (P7°) 2 - 2|pfr|p°] , (24) 

V + = \[Wi- Po° - Pr ) 2 - l^l 2 - (Pr°) 2 + 2|^b 7 °] . (25) 

Here, y is a dimensionless parameter defined as y = (pi — p ~ Pi) 2 1 s with = m z , and 
| Po| = (po ) 2 — m B c • The sum of the soft and hard sectors gives the total contribution of 
real corrections, i.e., T Real = r^ al + T h ^. 



With the real and virtual corrections, we then obtain the total decay width of Z boson 
to B c at the NLO accuracy of QCD 

^total r Born ^Virtual ^ Real + 0(at) . (26) 

In above expression, the decay width is UV and IR finite. In our calculation the FeynArts [l^| 
was used to generate the Feynman diagrams, the amplitudes were generated by the FeynCalc 



19(, and the LoopTools |20| was employed to calculate the Passarino-Veltman integrations. 



The numerical integrations of the phase space were performed by the MATHEMATICA. 



IV. NUMERICAL RESULTS 

To complete the numerical calculation, the following ordinarily accepted input parame- 
ters are taken into account: 

m b = 4.9 GeV, m c = 1.5 GeV, m z = 91 GeV, m w = 80 GeV, (27) 

Gp = 1.1660 x 10" 5 GeV- 2 , g 2 = 8Gp ™w = 0.4221, (28) 

v2 

in 



^ Bc (0) = ^=2 = 0.3616 GeVf . (29) 

V47T 



Here, R(0) is radial wave function at the origin of B c meson, which is estimated via potential 



model 



2lj and Gf is Fermi constant in weak interaction. The one loop result of strong 



coupling constant is taken into account in our calculation, i.e. 



47T 

= 7T, 9 — TTl ,,2 ^ • ( 30 ) 



QCD 



With the above preparation, one can readily obtain the decay width of Z° to B c meson 
in NLO accuracy of pQCD. In practice, the renormalization scale \i may run from 2m& to 
mz/2. At fj, = 2mb and then a s (2rrib) = 0.189 with Aqcd chosen to be 128 MeV, the LO 
and NLO decay widths of Z° — > B c + cb process are 

T LO (Z° ^B c + cb) = 72.31 keV (31) 

and 

T NLO (Z° -»■ B c + cb) = 78.45 keV, (32) 

respectively. And, at the scale \i = m z /2 and then a s {m z /2) = 0.140, the corresponding 
results are 

T LO (Z° ->■ B c + c6) = 39.43 keV (33) 

and 

T NLO (Z° ^ B c + cb) = 62.53 keV . (34) 
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Our LO result agrees with that existing in the literature Ref.|9j in case we take their 

3 

inputs, i.e. a s (m z /2) = 0.150 and ipB c (Q) = 0.369 GeV 5 . The above result indicates that 
at high energy scale, the NLO QCD correction to the decay width, or the B c production, is 
substantial. To see the scale dependence of the LO and NLO results, the decay width T(/i) 
and the ratio T(fi)/T(2mb) are shown in Figj6]for fi varying from 2rrib to mz/2. Calculation 
results indicate that after including the NLO QCD corrections, as expected the energy scale 
dependence of the decay width T(Z° — >■ B c + cb) is reduced evidently. 
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FIG. 6: The decay width r(/x)(left) and the ratio r(/x)/r(2m c ) (right) versus renormalization scale 
/i in Z° boson decays. 

V. SUMMARY AND CONCLUSIONS 

In this work we have calculated the inclusive decay width of Z° boson to B c meson at 
the NLO accuracy of perturbative QCD. Supposing that there will be copious Z° data in 
the future at the "Z-factory" , our results are helpful to the precise study of B c physics, and 
may also tell how well non-relativistic model works for B c system. 

Numerical results indicate that the NLO QCD correction slightly increases the LO result 
for the process Z° — > B c + c + b when ji is at the low energy scale of 2m?,, while it becomes 
huge, even comparable to the LO result, when \i runs to the scale of m z /2. We find that 
the energy scale dependence of the decay width is depressed, as it should be, when the 
next-to-leading order correction is taken into account, which means the uncertainties in the 
theoretical estimation are reduced. 
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